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Abstract 

We propose a new algorithm which works effectively in global updates in Monte 
Carlo study. We apply it to the quantum spin chain with next-nearest-neighbor 
interactions. We observe that Monte Carlo results are in excellent agreement with 
numerically exact ones obtained by the transfer matrix method. 
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The Monte Carlo method has been an indispensable tool for theoretial study in 
physics. An essential point of the method is to generate configurations according to 
their Boltzmann weights. For this purpose one usually starts with an initial config- 
uration and obtains a new configuration sequantially by some rule that guarantees 
the ergodicity and satisfies the condition of the detailed balance. Widely used rules 
are the heat bath method and the Metropolis algorithm. The heat bath mehtod is 
useful for systems which have a small number of candidates for updating because 
one has to calculate weights for all candidates to determine which one should make 
the new configuration. The Metropolis algorithm is suitable for complicated sys- 
tems since in this algorithm it is enough to know the weight just for one candidate 
to decide whether one should accept or reject that candidate for updating. The 
problem in the latter is that the efficiency of updating would be bad if the selection 
of the candidate is not appropriate. 

In previous papers [Q, Q for quantum Monte Carlo study using the Suzuki- Trotter 
formula we proposed re-strucuring method to improve the negative sign problem and 
applied it to one dimensional spin- 1/2 system with next-nearest-neighbor interac- 
tions. It turned out that this method is very useful for this model. We found 
however that it was difficult to make Monte Carlo simulations with large values of 
the Trotter number for technical reasons on the so-called global updates, which are 
necessary in the quantum Monte Carlo method in addition to the usual local up- 
dates. In the global updates one needs to change spin states on all sites along some 
lines connecting these sites, which usually depend the size of the system, while in 
the local updates limited part of the system irrelevant to whole size of the system 
should be updated. In the conventional apporach to the model, where the state on 
each site is represented spin up or spin down, conservation property of the quantum 
number (2;-component of the total spin J) is helpful to carry out global updates 
for large sizes of the system. In the re-structuring approach, on the contrary, we 
employ not eigenstates of al, where al denotes 2;-component of the Pauli matrix on 
site i, but eigenstates of a2i_i(T2i, that is 0jor — Ij. They are eigenstates of 

^2i~i + with eigenvalues 2, 0, 0, —2. Because of two Jz = states we need to con- 
sider ~ 2^" candidates for a global update in the Trotter direction with the Trotter 
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number n, which causes shortage of computer memory when n becomes large. Even 
for small values of n we found the Metropolis algorithm gives too poor efficiency to 
obtain results. In previous Monte Carlo work we therefore simulated systems with 
the Trotter number n < A only by employing the heat bath method with an approx- 
imation to discard AJz = 2 candidates. The agreement between Monte Carlo data 
and numerically exact calculations are satisfactory, but small discrepancy remains. 
In this paper we propose hybrid use of the heat bath method and the Metropolis 
algorithm which is applicable for large values of the Trotter number. We will see 
the discrepancy mentioned above is completely resolved. 

Now let us discribe the hybrid method. In this method we choose a candidate 
state applying the heat bath method repeatedly and use the Metropolis algorithm 
to determine whether we accept it or not. For concrete description we consider 
the quantum spin-1/2 system on an A^-sited chain with next-nearest-neigbor inter- 
actions. After the re-structuring the system reduces to an effective quantum spin 
system sitting on an A^/2-sited chain which has nearest-neigbor interactions only. 
Along the Trotter direction we have 2n sites with the Trotter number n. Each site 
on such two-dimensional lattice, where denotes the location in the space 
(Trotter) direction, is characterized by a spin state with S = 1, ®, or — 1. 
For one n-direction global updates with a fixed i we have to find a set of S(^ij-j, 
(j = 1,2, ...,2n), as a candidate. Let us denote S^ij) in the old configuration as aj 
and its candidate as jSj hereafter. First we consider the plaquette Qi in Fig. 1(a) 
only, where it is easy to pick up Pi and P2 according to the probability P^^\ 

p(i) = (71, A, 72, P2)/ E E ^(71, A, 72, P2), 

with 7j being the state of the nearest-neighbor site of the site {i,j) and 

W {(a, Cb, Cc Cd) being the Boltzmann weight for a plaquette depicted in Fig.l(b)[]. 

Next we choose (3^ in the plaquette Q2 according to the probability P^'^\ 

P(2) = iy(/32,72,/53,73)/EW^(/^2,72,C3,73). 
C3 

^We have 4'' = 256 values of W in this model. Actually, only 62 of them are non-zero. 
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Similarly we choose (j = 3,4,..., 2n—l) sequentially according to the probability 
p(-') which is, for odd j, 

pU) = W{^„P„j,+,,P,^,)/ W^(7„/?„7,+i,0+i) 

or, for even j, 

P(^-) = W^(/?„7„/3,+i,7,+i)/ E W^(/?„7„0+i,7,+i)- 

Cj+i 

Finally we compare the probability Pnew and Poia 

^ n E W^(/?2fc,72fc,C2fc+l,72fc+l) E W^(72fc+l,/?2fc+l,72fc+2,C2fc+2), 
k=l,n-l C2fe+i C2fe+2 

-Pold = W{a2n, l2n, "l, 7l) 

^ n E W^(«2fc,72fc,C2fc+l,72fc+l) E W^(72fc+l,a2fc+l,72fc+2,C2fc+2), 
fc=l,ri-l C2fc+1 C2fc+2 

and decide, following the Metropolis algorithm, whether we replace ai, 02, a2n 
by Pi, P2, ■■■,P2n or not. It is easy to see that this method satisfies the condition of 
the detailed balance. In Fig. 2 we plot Monte Carlo results on the R ratio obtained 
by this method, together with numerical exact calculations by the transfer matrix 
method. We perform simulations on an = 8 chain with the Trotter number up to 
16. We observe that the Monte Carlo data and the results from the transfer matrix 
method almost coincide for every value of n. It should be also noticed that the 
agreement is remarkable even for large values of the inverse temperature (3 because 
of much better efficiency of the updating than the efficiency in the previous heat 
bath method. 

To summarize we propose a hybrid method for global updates in the Monte Carlo 
algorithm. Through numerical work for a quantum spin system we show the method 
is quite powerful. Since the procedure in this method is quite general we expect it is 
applicable to many systems in the quantum Monte Carlo studies. We also expect the 
method is helpful in the classical statistical mechanics with the cluster algorithm^, 
where the global updates become important especially when the correlation length 
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is large. We would like to emphasize that the hybrid method enables us to choose 
appropriate candidates for the update at only a little cost of the computer time and 
the memory. 
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Figure Captions 
Figure 1 

(a) A schema of a part of the two-dimensional lattice with the Trotter number n — A 
to show how the candidate states are chosen for the n-direction global update in 
the hybrid method. The horizontal (vertical) arrow denotes the space (the Trotter) 
direction. 

(b) A plaquette with sites a, 6, c and d. The arrows are the same as in Fig. (a). 
Figure 2 

The R ratio, the ratio between (Z+ — Z_) and (Z+ + Z_) of the re-structured 
quantum spin system calculated by the Monte Carlo simulations with the hybrid 
method (mc) and by the transfer matrix method (tm) on an iV = 8 chain. Here 
Z+(Z_) denotes number of positively (negatively) weighted configurations generated 
in the simulation. Exact values of R are also plotted for comparison. Monte Carlo 
results in the figure with the Trotter number n = 2, 4, 8 (n = 16) is an average of 
a dozen of data calculated from 10000 (100000) configurations after 2000 (20000) 
sweeps for thermahzation. The errors are statistical ones only. 
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